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Future air traffic management systems are required to balance the conflicting objectives 
of maximizing safety and efficiency of traffic flows while minimizing the climate impact of 
aviation emissions and contrails. Integrating emission and climate models together with air 
traffic simulations improve the understanding of the complex interaction between the 
physical climate system, carbon and other greenhouse gas emissions and aviation activity. 
This paper integrates a national-level air traffic simulation and optimization capability with 
simple climate models and carbon cycle models, and climate metrics to assess the impact of 
aviation on climate. The capability can be used to make trade-offs between extra fuel cost 
and reduction in global surface temperature change. The parameters in the simulation can 
be used to evaluate the effect of various uncertainties in emission models and contrails and 
the impact of different decision horizons. Alternatively, the optimization results from the 
simulation can be used as inputs to other tools that monetize global climate impacts like the 
FAA’s Aviation Environmental Portfolio Management Tool for Impacts. 


I. Introduction 

The development and evaluation of concepts and technology to support future air traffic management systems 
require a hierarchy of models ranging from real-time simulations to extensive field evaluations. Air traffic 
simulation models such as Airspace Concept Evaluation System (ACES) 1 , Center Tracon Automation System 
(CTAS) 2 , Future Air traffic management Concept Evaluation Tool (FACET)’ and others are used to design air 
traffic systems balancing the conflicting objectives of maximizing safety, meeting future demands for airports and 
airspace and increase efficiency of traffic flows in the presence of uncertain weather. The impact of aviation 
emissions and contrails on climate imposes another constraint on the design of aircraft and aviation operations. The 
understanding of the complex interaction between the physical climate system, carbon and other greenhouse gas 
emissions and aviation activity can be improved by the development of integrated assessment models that include 
emission and climate models together with air traffic simulations. 

The impact of various greenhouse gas emissions on the climate depends on the amount of emission of each type 
of gas, the residual amount of gas in the atmosphere and the lifetime of the gas. The concentration of CCE in the 
atmosphere today, 390 parts per million (ppm), is far above the natural range of 172 to 300 ppm observed before the 
20th century and the increase is attributed to man-made activities 4 . CCE is the most important anthropogenic 
greenhouse gas and its impact on the climate as measured by radiative forcing is larger than all other emissions such 
as methane, nitrous oxide, ozone and water vapor. It is estimated that aviation is responsible for 2% of all 
anthropogenic CO? emissions. Oxides of nitrogen, commonly referred to as NOx, created by the high temperatures 
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in the aircraft engines affect the climate indirectly by affecting the distributions of ozone and methane. There is also 
concern about the impact of contrails and cirrus clouds due to air traffic. Contrails are clouds that are visible trails of 
water vapor made by the exhaust of aircraft engines. Persistent contrails reduce incoming solar radiation and 
outgoing thermal radiation in a way that accumulates heat'. The emission of water vapor from the aircraft is small 
relative to all contributions of water vapor to the atmosphere, but may be of concern at higher altitudes. 

The modeling of aircraft emissions and their interaction with each other to change the concentration levels of 
different gasses in the atmosphere and the resulting impact of the radiative forcing on the equilibrium of the Earth’s 
atmosphere is complex and requires the use of coupled atmosphere-ocean general circulation models together with 
three-dimensional models of carbon cycle and chemistry of other non-COj greenhouse gases. These models are 
computationally intensive and unsuitable for studies involving the generation of multiple scenarios. Simple emission 
and climate models, based on the input/output relations of linear systems, capture the fundamental emission to 
climate impact behavior by careful selection of key variables and their dynamics. 

The impact of various greenhouse gases depends on the total concentration, effect per unit change in atmospheric 
concentration and the spatial distribution of the gas. All these quantities are influenced by the lifetime of the gas and 
depend on the interval of assessment, which may vary from a few decades to a few centuries. Climate metrics 6 are 
aimed at providing a common scale to compare different greenhouse gases. If the metrics are to be used as a tool in 
developing and evaluating aviation operations, they should be transparent and easy to apply. Global Warming 
Potential (GWP) and Aggregate Global Temperature Potential (AGTP) are some of the commonly used metrics. 

This paper integrates a national-level air traffic simulation and optimization capability with simple climate 
models and carbon cycle models, and climate metrics to assess the impact of aviation on climate. Several authors 7 8 
have used linear climate models and metrics. However, the simulation and analysis capability described in this paper 
provides the ability to look at the impact of aviation operations on air traffic metrics such as airport arrival and 
departure rates, sector congestion and air traffic delays while providing the environmental impact metrics using 
linear models. 

Section II provides background for simulating air traffic using flight information and atmospheric and airspace 
data. Section III describes the aircraft fuel consumption and emission models. Section IV provides the model for 
diagnosing regions of airspace that are susceptible for persistent contrail formation. Section V explains the optimal 
trajectory generation with environmental constraints for cruising aircraft. Section VI models the climate response to 
aviation induced emission and contrails in terms global surface temperature change. Section VII describes metrics to 
measure the climate impact of aviation emissions and contrails with very different characteristics. Section VIII 
presents results based on the application of the trajectory optimization algorithm for calculating wind-optimal and 
contrails-avoidance routes and the climate response model that assesses climate response to air traffic along these 
routes. Conclusions and future work are described in Section IX. 

II. Air Traffic System Model 

The air traffic system model uses flight information and atmospheric and airspace data as inputs and simulates 
air traffic in the National Air Space (NAS). Air traffic is modeled using the Future Air Traffic Management 
Concepts Evaluation Tool (FACET). FACET is a flexible simulation environment for exploration, development, and 
evaluation of advanced Air Traffic Management (ATM) concepts. FACET models en route airspace operations over 
the contiguous United States. Core features include air traffic simulation, prediction, visualization, and playback of 
actual traffic data. Four-dimensional aircraft trajectories are modeled using spherical-earth equations. Aircraft are 
flown along flight plan or great circle routes as they climb, cruise, and descend according to performance models. 

Aircraft in FACET are modeled as parametric 3-dimensional trajectories using an inertial reference system. 
Values of the parameters vary depending on the aircraft. The equations of motion can be defined by the equations, 

x = V cos cp cos y + u(x,y), 
y = V simp cos y + v(x,y), 

E=(T -D)V/mg, 

h = V siny, (1) 

y = {L cos<j> - mg cos y) / mV, 

<p = Lsin0/mV cosy, 
m = -o(h,V,T), 
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where x is the downrange position, y is the cross track position, E is the energy height, h is the altitude, y is the flight 
path angle, <j> is the roll angle, q? is the heading angle, m is the aircraft mass, Th is the aircraft thrust, D is the drag 
and V is the airspeed. The fuel flow o is a function of thrust, altitude and airspeed. n(x, y) and v(x, y) are the x- and 
y-components of the wind. 

ITT. Aircraft Fuel Consumption and Emissions Model 

Eurocontrol’s Base of Aircraft Data Revision 3.6 (BADA) 9 provides the fuel consumption model for computing 
aircraft fuel burn. The following equation calculates fuel bum for aircraft 

f = t- SFC ■ Th , (2) 


where f is the fuel, t is elapsed time, Th is thrust, and SFC is the specific fuel consumption. 

The emission model is based on the Aviation Environmental Design Tool (AEDT) 1(1 developed by the Federal 
Aviation Administration (FAA). Six emissions are computed including C0 2 , H 2 0, SOx, CO, HC and NOx. 
Emissions of CO 2 , H 2 0 and SOx (modeled as S0 2 ) are directly proportional to fuel consumption". The emissions 
are computed by 


E C o 2 =EI COi -f = 3155- /, 

E Hi0 = EI Hi0 ' / = 1 237 • / , (3) 

E S o 2 - E SOi • / = 0- 8 • /. 

The terms E co ^ ,EI C0 , E HO ,EI HO and E so ,EI S0 are emissions and emission index of C0 2 , H 2 0 and 

S0 2 . The emissions are in grams, and fuel burn is in kilograms. Emissions of CO, HC and NOx are 
modeled through the use of the Boeing Fuel Flow Method 2 (BFFM2) 12 . The emissions are determined by 
aircraft engine type, altitude, speed, and fuel burn and the coefficients in the International Civil Aviation 
Organization (ICAO) emission data bank. The climate response model in the next section assesses the 
impact of aircraft emissions on climate change. 

IV. Persistent Contrails Formation Models 

The formation of contrails has been under investigation since 1919. 13 According to Appleman 14 , contrails are 
clouds of tiny ice particles that form when a mixture of warm engine exhaust gases and cold ambient air reaches 
saturation with respect to water. Contrails form in the regions of airspace that have ambient Relative Humidity with 
respect to Water (RHw) greater than a critical value, r contr 15 . Contrails can persist when the ambient air is 
supersaturated with respect to ice i.e. the environmental Relative Humidity with respect to Ice (RHi) is greater than 
100%. 2 In this study, the regions of airspace that have RHw greater than r contr and RHi greater than 100% are 
considered favorable to persistent contrails formation. More details for the contrails prediction model and the 
weather forecasting system, Rapid Update Cycle (RUC), are provided in a previous paper 16 . 

This model will be used to compute the regions in the U. S. airspace that facilitate persistent contrails formation. 
The blue, green and magenta polygons in Fig. 1 depict the areas at 37,000 feet above sea level in the U.S. national 
airspace where atmospheric conditions are favorable for persistent contrails formation at 6 a.m., 7a.m. and 8 a.m. 
EDT on May 24, 2007, respectively. The critical relative humidity and RHi values are computed using the RHw 
values, pressure, and temperature data obtained from RUC. Figure 1 shows that the location, shape and size of 
potential contrail regions vary with time. 

V. Aircraft Trajectories with Environmental Constraints 

Wind-optimal trajectories are fuel-efficient trajectories between any two city-pairs. They are also 
environmentally friendly since most emissions are directly proportional to fuel usage. However, if the trajectory 
goes through a region where the atmospheric conditions are favorable to persistent contrail formation, contrails 
create additional impact on the climate. The amount of contrail formation can be reduced either by flying around the 
potential contrail regions at the same cruise altitude (2D approach) or by changing the cruise altitude and going 
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around the potential contrail regions (3D approach). The amount of additional fuel consumed in flying a contrail 
reduction route instead of flying a wind optimal route represents increased cost to the airlines as well as a negative 
impact on the environment due to the resulting excess CO 2 emissions. Each additional kg of fuel produces 
approximately 3.15kg of CCE.The optimization module in the simulation generates trajectories that reduce travel 
time through contrail regions by imposing a penalty on travel through potential contrail areas 16 . 



Figure 1. Regions of airspace at 37,000 feet favorable for persistent contrails formation between 6 a.m. and 8 
a.m. on May 24, 2007. 

Three optimal trajectories from Chicago O’Hare airport (ORD) to Newark Liberty airport (EWR), corresponding 
to three different penalties on the contrail regions, are shown in Fig. 2 for flights with cruising altitude equal to 
34,000 feet. The cruising speed is assumed to be 400 nmi/hr (741 km/hr). The green arrows represent the wind 
directions, obtained from RUC, at 6 a.m. EDT on May 24, 2007. The arrow sizes are plotted in proportion to the 
wind magnitudes. The areas favorable to persistent contrails formation (RHi > 100%) are indicated by the blue 
polygons. C, is a weighting factor used to impose a penalty on a trajectory going through a potential contrail region. 
The wind-optimal trajectory is generated imposing no penalty, C r =0, for travel through contrail regions. Two 
optimal trajectories in additional to the wind-optimal route are also plotted in Fig. 2. The optimal route with C r =0.6 
completely avoids the contrail polygons near the departure airport. The optimal route with C r =0.2 only partially 
avoids the polygons but is shorter. Note that both routes travel through the blue polygon surrounding the destination 
where aircraft start to land. In this case, there is a tradeoff between flying a shorter route with more persistent 



Figure 2. Optimal trajectories at 34,000 feet from ORD to EWR with different penalties on contrail regions. 
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contrails formation versus flying a longer route with less persistent contrails formation. The performance of optimal 
trajectories is evaluated by investigating the total travel time and the time spent traveling through regions of 
persistent contrails formation. 


VI. Linear Emission and Climate Models 

This section models the climate response to aviation emission and contrails as outputs from a series of linear 
dynamic systems. The linear systems are generated against a background of concentration of various greenhouse 
gases resulting from past emissions from all sources. Linear emission models provide the incremental changes to the 
greenhouse gas concentrations resulting from the emission due to some or all aircraft operations. 

Several linear models have been proposed for representing the change in the concentration of CCL due to the 
emission of a unit amount of CCKThe carbon cycle models describe the changes to the CO 2 concentration due to the 
transport and absorption of CO 2 by the land mass and various ocean layers. A 19 th order model with states 
representing the interaction between different oceans and depths is described in the literature 17 . This paper uses a 3 rd 
order linear system to describe the CO 2 concentration dynamics and the model parameters will be described in 
subsequent sections. The concentration dynamics of other non-CCL greenhouse gases are described by first order 
linear systems. Radiative forcing (RF) associated with CO 2 and other gases can be adequately approximated by 
simple functions of concentrations 18 " 19 . 

Radiative Forcing due to different emissions affects the climate by changing the global average temperature. The 
temperature response/energy balance to RF can be modeled using either a first order linear modeL 0 or a second order 
linear model 21 . The time constants in the two-box ocean model correspond to the dynamics associated with the 
surface layers of the ocean and the thermal inertia associated with the deep ocean. A comparison of different linear 
models used to estimate the impact of aviation on climate suggests that multi-box ocean model gives closer 
approximation to ocean temperature response 7 . The modeling of emission to temperature change is illustrated in 
Figure 3. 
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Figure 3. Linear climate response models for aviation induced climate change. 

dx(t ) 


The response to a first order linear system 


dt 


■ = ax(t) + bu(t) where x(t) is the state variable and u(t) is the 


forcing function with a and b constants is given by the expression 


x{t)=e a( ' ‘ o) x(t 0 ) + J e a(t S) bu(C)d^. (4) 

? 0 
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Let A x(t) be the concentration change for a gas in the atmosphere at time t due to a disturbance u(t 0 ) applied at 
time, t 0 , resulting from a source of aviation emission. Since it is desirable to estimate the change in concentration, 
the initial value of x(t ) at t = t 0 can be assumed to be zero. This reduces the solution to 

t 

A x(t) = J e a(t ~ ■ ) bu(£,)d (5) 

The concentration change for a gas due to a unit impulse of gas emitted at t 0 is obtained by solving Eq. (5) and is 
equal to 


Ax(t) = be a( ' ‘ 0> . 


( 6 ) 


Similarly, the concentration change for a constant unit emission from t 0 to t is A x(t) = [1 -e aU ,0> ] . 

a 

The radiative forcing imposed by the emission can be assumed for the simplest model to be linearly proportional 
to the concentration change. The change in the radiative forcing is represented by 

t 

A F(t) = J* kb(t -£)Ax{L,)dL, = kAx(t) that drives a surface temperature change AT(t) . 

{ 0 

A first-order climate model’" proposed in the literature models the temperature change, A T(t) , as the following 

cf ^W__A7X0 + Af , (() (7) 

clt A 

Substituting a = - 1 / AC and b = l/C in Eq. (7), the change in the surface temperature at time t is given by 

t 

AT(t,t 0 ) = \/C f AF(t)d£. (8) 

f o 

The change in the surface temperature to an exponentially decaying radiative forcing of strength A, which is in the 
same form as shown in Eq. (8) and defined as A F(t) = Ae resulting from the emission applied impulsively 

at the initial time t 0 is 


AT(t,t 0 ) = 1/C f e H, - S)UC Ae~ (£ ~ to)la d^. (9) 

h 

Evaluating the integral, Eq. (9) reduces to 

A = Aar[e~ (t ~ ,o),a -e~ { '~‘ a)lr ]l C{a - r) (10) 

where the time scale of climate response T = AC . C is the heat capacity of the climate system and A is the climate 
sensitivity parameter. Climate response to aviation emission can be assessed through quantifying global surface 
temperature change modeled by cascading several first order systems. The following subsections introduce a 
recently developed climate assessment metric based on linear climate systems for modeling temperature change with 
higher order systems. 


VII. Climate Metrics 

Aviation emissions and contrails have very different characteristics and influence the climate either directly or 
by decomposition into other chemical compounds. Figure 4 shows the lifetime associated to contrails, ozone, and 
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CO 2 emission. The impact of certain gases depends both on the amount of emission and the location of the emission. 
In addition, the impact of the emissions varies with the choice of the time horizon. These variations make it 
necessary to develop a common yardstick to measure the impact of various gases. Several climate metrics have been 
developed to assess the impact of aviation emissions. These metrics have been used to develop strategies for limiting 
the impact on the environment 6 " 2 . 



years 


Figure 4 . Lifetime associated to contrails, ozone, and CO2 emission. 

A. Global Temperature Change Potential 

Global Temperature Change Potential (GTP) 23 is a climate assessment metric that adapts a linear system for 
modeling the global temperature response to aviation emissions and contrails. The definition of the absolute GTP is 
a convolution integral from t n = 0 to t=H, and has the following representation, 

H 

AGTP(H) = fR(H -DA F(£)d£, (11) 

0 

where R(H -£) is the impulse response function for the surface temperature change at time //due to a radiative 
forcing applied at t, . Note that surface temperature change AT(t,t 0 ) by the simplified climate model in Eq. (7-10) 

is equivalent to the AGTPfTZ) when setting t=H, t 0 = 0, R( H - '£ ) = ~ e ~ ( H " , , and AF(D = Ae " la . Two 

versions of AGTP are available in the literature. The pulse AGTP measures the change in the global temperature at 
a particular time, t, in the future due to an instantaneous disruption at t 0 . The sustained AGTP measures the global 
temperature change at time t due to disruptions constantly applied for a period between t and t 0 . 

AGTP is employed in this study for translating aviation induced CO 2 emission and persistent contrails into total 
effect on global warming. In general, higher order systems are adopted for modeling the impulse response function 
R(H - D and the radiative forcing AF(g). The formulations for AGTP due to CO 2 emission and contrails are 
provided in the following subsections. 

B. Pulse and Sustained AGTP for CO2 emission 

The response of global-mean surface temperature to a unit impulse of A F , using a second order model - is given 
by 


(id 

H C J 

where the parameters Cj and dj are given in Table 1 taken from the literature 6 . The first term and second term in the 
summation are associated with the response of the ocean mixed-layer and Deep Ocean, respectively, to a forcing. 
The radiative forcing for CO 2 emission is made of a steady state component and three exponentially decaying 
components, and is given by 
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( 13 ) 


A F c ° 2 (/) = A c ° 2 (a 0 + J a ie - (t - ,o)lai ) , 

i = 1 

where A c 2 =1.82xl0" 15 Wm' ] kg' 1 is the specific forcing due to COj and its value is taken from past studies 24 . The 
exponential component in parentheses is the decay of an impulse emission of CO? with time and the parameters are 
shown in Table l.The expression for AGTP given in Eq. (1 1) is evaluated after substituting Eq. (12-13). 

The pulse AGTP for 1 kg CO 2 emission for a time horizon H can then be found by analytical integration using 
Eq. (1 1-13), and is given by 


j - 1 


i=i 


AGTP C0 2 (H) = A c ° 2 ^ a o c j( 1 - e ~ H ' di ) + J — C 


'Me 

d. 


- H/a , ~ H/d i\ 

• -e J ) 


(14) 


This measures the surface temperature change due to CO? emitted at a particular time instant at t=0. 

The sustained AGTP for CO 2 models the temperature change due to CO 2 emitted constantly over a period of time 
between t=0 and t=H. For the same amount of CO 2 emissions, the sustained AGTP is formulated as 



(15) 


H-dj( \-e~ Hldj ) +2] 


OjOfj 


OC: 


_ /1 — H / cc : \ j /1 — Hid: x 

aAl-e ‘)-dA 1-e 1 ) 


Table 2 shows the pulse and sustained AGTP values per kilogram of CO 2 emission for three time horizons. The 

absolute global temperature potentials are expressed in terms of the specific forcing, A c 2 , since different values 
can be selected from the literature. 


C. Pulse and Sustained AGTP for Persistent Contrails 

The surface temperature response for contrails is modeled similarly. An impulse function best characterizes 
contrails radiative forcing since contrails are short-lived, usually last for several hours, in the atmosphere. The pulse 
AGTP for contrails formation is simply taken as the impulse response. For a unit of contrails induced radiative 
forcing, the pulse AGTP is represented by 


H 

AGTP CoDtnil \ H ) = f R(H - r C)d( r C-(Y)d r C = R(H). 

0 

The sustained version for AQ , pp Contlalls obtained by integrating R(t) over the period of H, 


1 

H 


H 2 

0 j. 1 


-Hid, 


) 


(16) 


(17) 


Equation (16-17) computes surface temperature change to a unit of energy induced to the atmosphere by 
contrails. The net radiative forcing for contrails includes the longwave RF LW and the shortwave RF SW radiative 
forcing and is defined as RF nels = RF LW + RF SW . It is measured in terms of unit of power (Watts) per unit area of 
contrails (m 2 ). Typical values for RF nets have a range between 10 Wm" 2 and 30 Wm' 2 taken from Meerkotter 5 and 
Haywood 25 . Due to the nature of contrail formation, it is argued that it is better to represent contrail radiative forcing 
in terms of unit distance flown by the aircraft (watts/km). The amount of energy, EF, induced to the atmosphere for 
a unit length of contrail over its lifetime is defined as 26 
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( 18 ) 


EF= f RF nets (Z)WM)dt , 

lifetime 


where W c is contrail width (m). Suppose RF net = 10Wm' 2 , contrail width )T=1000m and contrails lifetime is 
10000s, the energy EF for a km contrail equals 100 Gigajoules (GJ). The global surface temperature change, 
AGTP 0 mm ' h (H), is then computed by multiplying either Eq. (16) or Eq. (17) by the energy after it is normalized 
by the surface area of the Earth and total seconds in a year. Assuming iTE^lOOGJ per kilometer of contrails, the 
pulse and sustained AGTP values are computed for 3 time horizons using equations 12,16 and 17. The AGTP values 
are shown in Table 2. Note that they can be linearly scaled when other RF nets values are chosen. 

Table 1. Parameter values for the radiative forcing function due to CO 2 emission and impulse response 
functions in AGTP. 


Index 

0 

1 

2 

3 

a i (unitless) 

0.217 

0.259 

0.338 

0.186 

a j (Years) 


172.9 

18.51 

1.186 

c j (K/WnT 2 ) 


0.631 

0.429 


dj (years) 


8.4 

409.5 



T able 2. Absolute global temperature potential values for CO 2 emission and contrail s. 


Horizon (years) 

H=25 

H=50 

H=100 

Pulse AGTP C ° 2 (H) 
(K/kg) 

A C ° 2 x 0.3686 

A c ° 2 x 0.3176 

A c ° 2 x 0.2788 

Sustained AGTP C ° 2 (H ) 
(K/kg) 

A C ° 2 x 0.3024 

A c ° 2 x 0.3223 

A c ° 2 x 0.3079 

Pulse AGTP Conm ' h (H) 
(K/km) 

2.9936e-14 

6.9778e-15 

5.1045e-15 

Sustained AGTP Contm ' s (H) 
(K/km) 

1.5522e-13 

8.4377e-14 

4.5003e-14 


CO 2 is the most pervasive of all anthropogenic emissions. The impact of the introduction of additional amount of 
CO 2 on climate is better understood than the impact of all other greenhouse gases and contrails. In the development 
of climate metrics, the impact of other emissions is normalized with respect to the impact of CXX The atmosphere 
quickly absorbs additional CO 2 and the CO 2 radiative forcing is due to the globally distributed energy resulting from 
the uniform increase in CO 2 concentration. Contrails occur at different regions of the earth and add non-uniform 
sources of energy to the atmosphere. The efficacy factor is used to differentiate the way radiative forcing from CO 2 
and contrails affect the climate. Specifically, it is defined as the ratio between the global temperature increase for a 
local energy input relative to that for a CCVequivalent globally distributed energy input 26 . The efficacy factor for 
annual mean contrail cover is estimated to be about 0.6 27 and treated as a parameter in this paper ranging between 
0.6 andl.0. 


VIII. Results 

This section presents results based on using the simulation and analysis capability described in the previous 
sections. The optimal trajectory algorithm is applied to calculate an aircraft trajectory in the presence of winds that 
minimizes fuel burn and avoids regions of airspace that facilitate persistent contrails formation. The potential 
impacts to the climate for the U.S. domestic flights between 12 origin-destination pairs are assessed in terms of 
mean surface temperature change, i.e., AGTP due to aircraft emissions and persistent contrails formation. The same 
city-pairs were used by the Federal Aviation Administration to assess the impact of implementation of Reduced 
Vertical Separation Minima (RVSM) on aircraft-related fuel bum and emissions 28 . This study adapts the standard in 
RVSM and assumes that the cruising altitudes are between 29,000 and 41,000 feet. Eastbound aircraft fly odd 
thousands of feet while westbound traffic fly even thousands of feet. Figure 5 shows the wind-optimal trajectories 
for the eastbound flights at 37,000 feet at 6 a.m. Eastern Daylight Time (EDT) on May 24, 2007. Blue polygons 
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depict the areas favorable to persistent contrails formation. The trajectory computations are done using traffic and 
atmospheric data in the continental United States for May 4, May 24, and May 27 in 2007. The data for wind speed 
and direction are obtained from RUC. 
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Figure 5. The wind-optimal trajectories for the eastbound flights for 12 city pairs at 37,000 feet, 6 a.m. EDT 
on May 24, 2007. 


Persistent Contrails Formation in May, 2007 

Great-circle trajectories between the 12 city pairs are simulated to assess the severity of potential contrail 
formation induced by the air traffic throughout May in 2007. The results, shown in Fig. 6, is used to identify the 
days in May, 2007 that have low, medium and high potential contrail formation for the 12 city pairs. Then, climate 
impact of air traffic is investigated in more detail for these days. Persistent contrail-favorable regions over the 
continental United States are computed 16 for May 2007 using atmospheric data from the RUC. A flight along the 
great-circle route at the assigned flight-plan altitude between the 12 city pairs is simulated for each hour thorough 
May 2007 since atmospheric data are updated hourly. The total time that aircraft fly inside contrail-favorable regions 
is recorded each day. Figure 6 shows the average contrail formation time for a group of 24 flights traveling between 
the 12 city pairs. May 24 has low (55 minutes) contrail formation. May 4 has medium contrails formation, and May 
27 has the highest (493 minutes) contrail formation. 
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Figure 6. Average persistent contrail formation time for each day in May 2007. 
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Trade-off between AGTP and Additional Fuel Burn 

This section analyzes climate impact of the wind-optimal and contrails-avoidance trajectories for 12 origin- 
destination pairs on low, medium, and high contrails formation day. The optimal aircraft trajectories are generated at 
the beginning of each hour for a period of 24 hours starting from 6 a.m. EDT on May 24, 4, and 27 in 2007, 
respectively, using hourly updated weather data from RUC. Six flight levels are considered for each direction of air 
traffic for each city pair. A group of 21 optimal aircraft trajectories are calculated for each flight level by increasing 
the value of C r from 0 to 2 with increments equal to 0.1. The cost coefficient of time is chosen as C, = 20 for each 
case. The cruising speed is assumed to be 420 nmi/hr. The fuel consumption for each aircraft trajectory is 
calculated using Base of Aircraft Data 9 (BADA) formulas by assuming that the aircraft are short to medium range 
jet airliners with medium weight. CO 2 emission is then obtained from fuel burn using the emission index. The 
persistent contrails formation time associated with each trajectory is also recorded. The details of aircraft trajectory 
generation are described in an earlier paper 16 . 

A total of fifty bins are defined such that the aircraft trajectories can be categorized based on their additional fuel 
consumption. The first bin contains the wind-optimal trajectory, which is the baseline for fuel use comparison and 
corresponds to trajectories that require zero % of additional fuel consumption. The second bin contains aircraft 
trajectories that consume less than 1% additional fuel, the third bin contains those consuming less than 2 %, and etc. 
The fiftieth bin has trajectories that burn more than 49% of fuel. In each bin, the optimal trajectory that has least 
amount of persistent contrails formation time is selected to represent the bin. Note that there are fifty bins for each 
group of trajectories and six groups for each direction of air traffic every hour. 

The surface temperature change due to CO 2 emission is calculated for the trajectories in each bin. Figure 7 shows 
the AGTP values and fuel burn for flights between the 12 city pairs on May 27, 2007. The average fuel bum for 
each set of 24 flights along the wind-optimal trajectories is about 130 tons that is converted to three AGTP values 
for time horizons of 25 years, 50 years, and 100 years. The climate response (AGTP) to CO 2 emission diminishes as 
the time horizon increases and increases linearly with the fuel burn for each horizon. 
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Figure 7. Surface temperature change due to fuel burn for May 27, 2007. 


There is significant uncertainties 29 " 31 with the radiative forcing associated with contrails and the equivalent 
energy forcing. The surface temperature changes due to contrails formation associated with three different radiative 
forcing values, which are 30 Wm" 2 , 10 Wm" 2 , and 3.3 Wm' 2 are plotted in magenta, blue, and green in Figure 8. The 
dash-dot and solid curve show results for flights with and without altitude optimization, respectively. The error bars 
indicate the variation of AGTP due to contrails as the efficacy factor is varied from 1 to 0.6. The figure shows when 
altitude is optimized, about 6 tons increase in total fuel consumption can reduce the surface temperature change due 
to contrails to OK from 5.7xlO" 10 K, 1.9xlO" 10 K, 0.63xl0" 10 K for various radiative forcing values when time horizon 
is 25 years. Without altitude optimization, the reduction in AGTP values due to contrails is more gradual with 
increase in total fuel consumption. Allowing further increase in fuel consumption does not result in proportionate 
reduction in contrail induced AGTP. The AGTP values are smaller but the trends remain the same for these curves 
when a longer time horizon is considered. 
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Figure 8. Surface temperature change due to contrails formation for May 27, 2007. 


Figure 9 shows the total surface temperature response to both the CO 2 emissions and contrails. Without 
optimizing the altitudes, the set of trajectories that burns 136 tons of fuel acquire the minimum total surface 
temperature change, 6.79xlO' 10 K, for 25 years time horizon when contrails radiative forcing value is 30Wm' 2 . The 
sets of trajectories that bum 135 tons and 133 tons have the smallest AGTP values, 4.17xlO" 10 K and 3.28xlO" 10 K, 
when the radiative forcing are lOWm" 2 and 3.3Wm" 2 . Burning extra 4.6%, 3.9%, and 2.3% fuel reduce the AGTPs 
by 20%, 11%, and 3.2%, for the 3 radiative forcing values for 25 years time horizon. Aircraft trajectories selected 
for minimum temperature change are closer to wind-optimal when contrails radiative forcing is smaller. Because 
wind-optimal trajectories minimize CO 2 emission that has a life time much longer than contrails, similar trends are 
observed when a longer time horizon is considered. The minimum AGTP values are 2.89xlO" 10 K, 2.88xlO' 10 K, and 
2.87xlO" 10 K for the three contrails radiative forcing when attitude is optimized. The associated fuel burn values are 
about 137 tons, 136 tons and 136 tons. In this case, burning extra 5.4%, 4.6%, and 4.6% fuel reduce the AGTPs by 
66%, 38%, and 15%, respectively. When 100 years time horizon is considered, 4.6%, 4.6%, and 3.1% extra fuel 
burn reduces the AGTPs by 29%, 10%, and 1.8%, respectively. The temperature changes are almost entirely due to 
CO 2 emission since these trajectories almost completely avoid contrail favorable areas. Aircraft, with optimized 
altitude, are also flying closer to wind-optimal for minimum temperature change when contrails radiative forcing is 
smaller or horizon is longer for the same aforementioned reasons. 

H=25yrs H=100yrs 
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Figure 9. Total surface temperature change due to fuel burn and contrail formation for May 27, 2007. 


Figure 10 and Figure 1 1 show the surface temperature change for a medium and low contrails day, respectively. 
The figures show that the reductions of surface temperature change are smaller for these days because contrails are 
less severe. This is particularly true when time horizon is longer where AGTP are mostly contributed by CO 2 
emissions only. The curves for May 24 are almost flat since contrails are very low. Flying wind-optimal should 
minimize climate change and the fuel cost. These curves can be used to develop optimal fleet allocation and 
scheduling strategies to minimize climate change for a given amount of additional fuel burn. 
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Figure 10. Total surface temperature change due to fuel burn contrails formation for May 4, 2007. 
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Figure 11. Total surface temperature change due to fuel burn contrails formation for May 24, 2007. 


IX. Conclusion 

This paper described a simulation capability that can be used to evaluate technology and concepts to mitigate the 
impact of aviation emissions on climate. The capability brings together metrics useful in aviation operations with 
those used in climate studies. The capability can be used to make trade-offs between extra fuel cost and reduction in 
climate impact. When altitude is optimized during a severe contrails day, burning extra 5.4%, 4.6%, and 4.6% fuel 
reduces the global surface temperature change by 66%, 38%, and 15%, respectively, for contrails radiative forcing 
of 30Wm' 2 , 10Wm' 2 , and 3.3Wm' 2 for a 25-year time horizon. When 100-year time horizon is considered, 4.6%, 
4.6%, and 3.1% extra fuel reduces the absolute global temperature potentials by 29%, 10%, and 1.8%, respectively. 
Aircraft trajectories selected for minimum temperature changes are closer to wind-optimal because they minimize 
CO 2 emission that has a life time much longer than contrails. This condition also prevails when contrails radiative 
forcing is smaller or a longer time horizon is considered. Similar trends are observed during the medium and low 
contrails days in which the reductions of surface temperature change are smaller because contrails formations are 
less severe. There is considerable uncertainty in our understanding of the radiative forcing associated with 
emissions and contrails. The parameters in the simulation can be used to evaluate the effect of various uncertainties 
in emission models and contrails. It can also be used to evaluate the impact of different decision horizons. 
Alternatively, the optimization results from the simulation can be used as inputs to other tools that monetize global 
climate impacts like the FAA’s Aviation Environmental Portfolio Management Tool for Impacts. 
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